Numerical computation of one-photon mazer resonances for arbitrary field modes 
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We present a novel approach for solving numerically one-dimensional scattering problems and 
apply it for computing the emission probability of an ultracold atom interacting with an arbitrary 
field mode of a high-Q cavity. Our method is efficient, stable and succeeds when other numerical 
integration methods fail. It also replaces and improves advantageously the WKB approximation. 
The cases of sinusoidal, sech 2 and Gaussian field modes are studied and compared. Divergences 
with previous works, where WKB was used, are found. 
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I. INTRODUCTION 

Some recent works have been devoted to the interac- 
tion of ultracold atoms with microwave cavities 
These studies treated the interaction between an incident 
atom in an excited state and a cavity field containing n 
photons, taking the quantum mechanical CM motion of 
the atom into account. This interaction leads to a new 
kind of induced emission intimately associated with the 
quantization of the CM motion, named the mazer action 

After interaction with the cavity, the atom can be 
found transmitted in the excited state or in the lower 
state, or reflected as well. It has been shown that 
these event probability amplitudes (denoted respectively 
T a , ni T bn+ i, R a ,n and Rb, n +i) are given by 



Pa,n 2 t—n) ^b,n+l 

Ra,n = \{r+n + r- n ) Rb. 
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where t± n and r± n are, respectively, the transmission and 
reflection amplitudes of the elementary scattering process 
of the atom incident upon the potential 



V=(x) = ±hg\/n+ 1 u(x) 



(2) 



where x is the atom traveling direction, g is the atom - 
field coupling strength and u(x) is the cavity field mode 
function. 

The induced emission probability of the atom interac- 
ting with the cavity field is then given by 



Pern — |Tfc.„+i| 2 + |i?i,, n +l| 



(3) 



All these probabilities are strongly dependent on the 
cavity mode profile. Their calculation needs to solve the 
one-dimensional time-independent Schrodinger equation. 
Loftier et al. || have calculated them as a function of 
the interaction length n n L (L is the cavity length and 
n n = K\/n + 1 with k = y/2mg/h, and m the atomic 



mass) for various cavity mode functions: mesa, sech 
and sinusoidal modes. For the two first modes, analyti- 
cal results of the probability P ern (n n L) have been given. 
For the sinusoidal modes, detailed Wentzel-Kramers- 
Brillouin solutions have been presented. Nevertheless, 
Retamal et al. || have shown that the WKB approxima- 
tion may lead to inaccurate predictions for actual inter- 
action and cavity parameters. 

The consideration of actual interaction and cavity pa- 
rameters is a difficult numerical task as they can corre- 
spond to very large values of K n L. When no analytical 
solutions of the Schrodinger equation are available, usual 
numerical integration methods do not converge rapidly 
or do not converge at all. 

In this paper, we present a novel approach for solving 
the one-dimensional time-independent Schrodinger equa- 
tion. Our method is described in Sec. |n]. It may be ap- 
plied for any potential function. It is efficient, stable and 
succeeds when other numerical integration methods fail. 
It may replace and improve advantageously the WKB 
approach. Details of the implementation of our method 



and validity tests are presented in Sec. [II. 

We apply our method for calculating some new 
Pemi^nL) curves. In particular, the Gaussian potential 
has been considered, thinking in open cavities in the mi- 
crowave or optical field regime. These results are pre- 
sented in Sec. 
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II. DESCRIPTION OF THE METHOD 

The one-dimensional time-independent Schrodinger 
equation can be written in atomic units in the form 



2dx 2 



+ (E-V{x)) <p(x) = 



(4) 



We assume here that the potential V(x) has non 
zero values only in a given region of the x axis, say 
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[x a ,Xb\. This region is divided into J grid points, de- 
noted Xt,...,xj such that x a — x\ < X2 < ■ ■ ■ < x,j-i < 
Xj = Xb- Let us note Ij (0 < j < J) the region [xj,Xj + {\, 
I ]-oq,Xi], and Ij [xj,+oa[. 

The potential V(x) is approximated in each region 
Ij by a straight line connecting V(xj) and V(xj+x). 
The approximated potential is noted V a p prox (x). Then 
Schrodinger equation takes the simple form in Ij (0 < 
3<J): 



(5) 



with ao = a,/ = 2E, bo — bj — 0, and for < j < J 



a =2 (e- x o+xV j -x i v j+l \ 



6 , = - 2 Srt^t with v i = v(x i ),i = j,j + i 



(6) 



The most general solutions of Eq. (g) are given by 

$ j {x) = C j f+{x)+D j fr( x ) (7) 

where Cj and Dj are two complex constants and fj(x) 

and fj~(x) 2 functions depending on the a,j and bj values. 
These functions are given in Table ||. 

As shown in Table |, the functions ft(x) and f7(x) 
depend on the sign of cij + bjx. To avoid a change of this 
sign in a given region Ij , the set x\ , . . . , xj must contain 
the roots of the equation V(x) = E. Thus, the condition 
a,j + bjx > (resp. a 3 + bjx < 0) is equivalent to that 
E > V(x) (resp. E < V(x)) for x £ Ij. 

The complex constants Cj and Dj in Eq. (Q) are deter- 
mined by use of the wavefunction asymptotic behaviour 
knowledge ((C 0} D ) or (Cj,Dj) are supposed to be 
known) and of the set of conditions imposing the con- 
tinuity of the wavefunction and its derivative along the 
x axis : 



(j>j{Xj) = 4>j-l(Xj) , , 

= #-ifo) ( " <J - J >- 



If (Co, Dq) are given, Eq. (g) may be written as 



(8) 



%)= 1 » J > (.!» 



with 



/j " fj Hi 
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and 



gt(x) = ^,i = j,j-l. (11) 

If (Cj,Dj) are known, Eq. (0) should be written 



D; 



D , 



for J > j > 0, (12) 



with 



fj9 S -x - fj if fjjj-x - fj-iQj 
fj-i9j - fj 9j-i fj-i9j - fj 9j-i 



(13) 



The knowledge of 4>j{x) in each region Ij defines an 
approximated wavefunction of the Schrodinger equation 
(Q). The grid point number J fixes the accuracy of the 
method. The higher J is, the more accurate is the ap- 
proximated potential V approx (x) and the better is the 
approximated wavefunction. The actual number J to 
be considered for a given accuracy depends on various 
parameters (potential form, ener gy E , size of the region 
[x a , Xb])- We show in Sec. [II and IV that a few hundreds 
are typical numbers. 

The wavefunctions yielded by our method are not di- 
vergent at the classical turning points (where E — V(x)), 
on the contrary of the WKB method. Although some 
Bessel functions of Table B are divergent at these points 
(as Zj(x) = cij + bjx = 0), all the functions f^ix) and 
their derivative gf(x) admit finite limits there (see Ta- 
ble ||. 

The transmission and reflection complex amplitudes 
(denoted respectively t and r) of a particle incident upon 
a potential V(x) may be calculated as well (and con- 
sequently the induced emission probability (||)). If we 
consider the outgoing wavefunction (f>j(x) = e lkx with 
k = V2E, we can calculate the corresponding incoming 
wavefunction 4>q(x) — Co cos(fca;) + Dq sin(fca;) by use of 
relations (113) . We then have 



Co - iDn 



Cp + iP 
Co - iD 



(14) 



III. METHOD IMPLEMENTATION 

We have implemented the method described above 
on a PC Pentium based computer. Bessel functions 
of Table Q have been coded using algorithms given by 
Zhang and Jin 0. The evaluation of the Mathemat- 
ica ® Bessel functions has been discarded due to poor 
performances (thousand times slower than direct compu- 
tation). For usual cavity parameters, Co and Dq coeffi- 
cients of Eq. (|lj) can rapidly become very large. To avoid 
computation overflow errors, we have adopted a custom 
number representation. Our internal range of values was 
9.9 E ± (~9 E 18) with 15 significants. 

As a first test of our method, we checked its ability 
to converge when the grid point number J increases. 
We have calculated the induced emission probability P em 
of an atom interacting with a sech 2 mode profile cavity 
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(u(x) = sech 2 (x/L)) as a function of the number J and 
for a fixed interaction length n n L. Fig. [l] presents the 
curve obtained in the case of k/ K n — 0.01 and K n L = 10. 
This curve shows very clearly that P em converges to a 
given value as J increases. Also, in that case, a J num- 
ber of 200 is sufficient for predictions to the accuracy of 
the graphics. For this calculation, the region [x a ,Xb] was 
limited to 16 times the cavity length L (~ 8 times the 
FWHM of the potential function). 

In order to improve the convergence of our method 
with the number J, we renormalized the approximated 
potential so that the area under it remained constant 
whatever value J had. Instead of using Vi — V(xi) in 
Eq. (^), we have used V% = aV(xi) with a equal to the ra- 
tio of the area under V(x) over the area under V a p prox (x). 

The study of the sech 2 mode profile is a good test for a 
numerical method as there exists an analytical expression 
for P em (n n L). In this case we have (see Loffler et al. ||]) 

= r [1/2 - i{kL + g )] ■ r [1/2 - z{kL - g*)] 

r {-ikL] • r [i - ikL] 

(15a) 

± _ T[ikL] ■ T[l-ikL] ± 

r " - r[i/2 + i£] • r [1/2 

With £ = y/±(K n L)*-l/4. 

Fig. H shows two P em (n n L) curves calculated for 
k/n n = 0.01 and k/n n — 0.1 in the case of the sech 2 
mode profile. On this figure, the solid lines represent the 
analytical results deduced from formulas ([!]) , (||) and ( |l5| ) 
while the dotted ones represent those obtained using our 
method (the grid point number J was fixed to 200). The 
agreement between these results is very good. 

IV. NEW RESULTS 
A. Fundamental sinusoidal mode 

The fundamental sinusoidal mode profile 

u (x) = { sin ( wx / L ) for 0<x<L 
v ' 1 elsewhere v ' 

has been studied by Loffler et al. [|| and Retamal 
et al. . Their conclusions are not in agreement with re- 
gard to the behaviour of the induced emission probability 
for ultracold atoms (k/n n — 0.01) interacting with large 
cavities (n n L of the order of 10 5 ). Such cavity parame- 
ters are of a great interest as they correspond to realis- 
tic values for Rydberg atoms interacting with microwave 
cavities in recent experiments performed by the Ecole 
Normale Superieure Group || (see discussion about the 
orders of magnitude in ||). Retamal et al. || have pre- 
dicted well resolved resonances in the curve P em (K n L) for 
the parameters cited above, whereas the curve predicted 



on the basis of the results of Loffler et al. looks very 
different. 

We have calculated the P em {K n L) curve for an atom 
interacting with this mode profile by use of our method. 
Fig. U presents this curve for n n L comprised between 
100000 and 100010 (fc/«„ = 0.01 and J = 100). Our 
curve shows the same well resolved resonances as those 
predicted by Retamal et al. || . So we confirm their pre- 
dictions that the WKB approximation cannot be used 
when considering ultracold atoms with k/n n — 0.01 and 
realistic interaction lengths of the order of 10 5 . Our 
calculations confirm also that these resonances have got 
smeared at K n L = 10 5 if we consider warmer atoms cha- 
racterized by k/n n = 0.1. For this last value of k/n n 
the WKB approximation still holds and P em is then in- 
variably equal to 1/2. It is good to remark here the 
importance of still having resolved resonances in P em for 
actual interaction and cavity parameters, as this result 
will permit the possibility of testing, in principle, the 
mazer effect with present cavities. The construction of a 
re-entrant cavity, as it was proposed in Ref. || , would not 
be necessary for testing the mazer effect, but for improv- 
ing the resolution of the resonances by approximating a 
mesa function mode profile. 

B. First excited sinusoidal mode 

The first excited sinusoidal mode profile 

u(x) = { sin ( 2 W L ) for < x < L (17) 
^ ' | elsewhere ^ ' 

has been studied by Loffler et al. || . These authors have 
derived an expression for P em on the basis of a detailed 
WKB calculation. They have shown that, for large in- 
teraction lengths (without being very explicit about the 
word "large"), the behaviour of the induced emission 
probability is described by P em (n n L) — sin 2 (A„) with 

1 W2 

A n = n n L ■ — / v(fc/K„) 2 + cos(x)dx (18) 
tt Jo 

which can be written A„ = K n L ■ Const, if k/ n n is a fixed 
ratio. 

In Fig. [|, we present two curves of P em (n n L) calcu- 
lated for K n L comprised between 100000 and 100020 (for 
k/n n = 0.01 and k/n n = 0.1). The first curve does 
not exhibit a square sine dependence over the interac- 
tion length, whereas the second one does. This indicates 
that K n L = 10 5 cannot be considered as "large" in the 
sense of Loffler et al. § for k/n n = 0.01. For kj n n = 0.1, 
A„ = 2TTK n L/T, with the period T equals to ^16.3 (ac- 
cording to Eq. (p^l)). This period is well reproduced by 
our calculations (see curve (b) on Fig. |^) . 

It is a remarkable result of Loffler et al. ||] the possibil- 
ity of building state-changing and state-preserving mir- 
rors for atoms by modifying the length of the cavity using 
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this first excited sinusoidal mode. From Fig. |4| we con- 
clude that resonances in P em for the first excited mode 
are even narrower than those predicted by a square sine 
function sin 2 (A„) . This is a very convenient result as we 
are considering actual interaction and cavity parameters. 

C. Gaussian mode 

Up to-date, the Gaussian cavity mode profile has not 
been studied exactly in the quantum theory of the mazer. 
Loffler et al. [|| have argued that the sech 2 mode profile 
could be used as a good approximation of the Gaussian 
one. 

To verify this assumption, we have calculated various 
Pem(K n L) curves for the profile 

u(x)=e~^ (19) 

The parameter a was fixed to ^J^JitL in order to adopt 
the same normalization factor for the two profiles (iden- 
tical area under the modes). 

Fig. |B| shows our results for k/ n n — 0.1 in the range 
n n L = to K n L = 20. Qualitatively both profiles ex- 
hibit the same behaviour : the resonances in the curves 
get smeared with increasing values of n n L. But this phe- 
nomenon is not so marked in the case of the Gaussian pro- 
file. Resonances still exist for longer interaction lengths. 
This is not a surprising result as the Gaussian profile is 
growing more abruptly than the sech 2 one. Thus it is 
in some sense "closer" to the mesa mode, which exhibits 
resonances at infinity. 

We have also considered the case k/K n = 0.01. Our 
calculations have shown that 90% damped oscillations 
are still present in the P em (n n L) curve for interaction 
lengths approximately 3 times larger in comparison with 
the sech 2 mode case. For these calculations, the region 
[x a ,Xb] was limited to 16 times the cavity length L and 
the grid point number J was fixed to 300. 

As it was pointed out in Ref. || the mazer effect is 
not restricted to the microwave domain and it might be 
tested more efficiently in the optical domain. Note that 
typical optical cavities have a Gaussian mode profile and 
it is possible to consider large coupling constants, fact 
that will help for testing the mazer effect. We want to 
call attention to the fact that Hood et al., in Ref. 0, 
presented the first experimental result for which the in- 
teraction energy Kg is greater than the atomic kinetic 
energy (k/n n < 1). 

V. SUMMARY 

We have developed a new method for solving one- 
dimensional scattering problems that may be applied ad- 
vantageously instead of the WKB approach. This has 
enabled us to calculate efficiently the induced emission 



probability P em of an atom interacting with a high-Q cav- 
ity for various mode profiles. Two sinusoidal modes have 
been considered. For these cases, we have been able to as- 
sert that the WKB approximation cannot be used in the 
computation of P em for ultracold atoms (kj ' n n = 0.01) in- 
teracting with actual cavities characterized by K n L of the 
order of 10 5 . Significant and convenient different physi- 
cal predictions are found, if we compare our results with 
previous works Q . 

The Gaussian mode profile has also been considered 
and we have shown that, although it exhibits a similar 
behaviour in comparison with the sech 2 mode profile, the 
resonances in the P em (K n L) curves exist for significantly 
larger values of n n L. The Gaussian mode is relevant 
when considering open cavities in the microwave or opti- 
cal domains. 

The presented numerical method will be helpful for 
computing the induced emission probability in the case of 
the recently studied two-photon mazer || , when consid- 
ering field mode profiles different from the mesa function. 
This and other applications of this numerical method will 
be presented elsewhere. 
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TABLE I. fj~(x) and f7{x) functions defining the general 
solutions (Q). In this table, J\/z(x) and Y\/$(x) denote respec- 
tively the first and second kind Bessel functions of order 1/3, 
Ii/ 3 (x) and Ki/z(x) the first and second kind modified Bessel 
functions of order 1/3, kj = sj2(E - Vjj, pj = \Jl(y 3 - E), 
and Zj(x) = a.j + bjx. 



b j= 


a>j = 


//(*) = 1 
/r (as) = x 


a,j > 


fj~(x) = cos(kjx) 
fy{x) = sin(fcja;) 


CLj < 


/+(*) = e-'i' 
fr(x) = eV" 




Zj(x) > 


tf(x) = y/\z j (x)\Ji 
ff(x) = ^\z j (x)\Y i 






Zj(x) < 


f;(x) = ^\ Zj (x)\n( 

fr{x) = y/\zj(x)\Ki 

J 3 




) 



TABLE II. Behaviour of the fj(x) functions and their 
derivative gf{x) at the classical turning points when bj 7^ 0. 



E > K(a;) 


region 


E < V(x) 


region 







/ + (*)-» 





s~0«0 -> 


--\3bA 1/3 T(l/3) 
n 

bj 


g + (x) - 

9~{x) -> 


iisfc.-l^rci/s) 


|36 J j 1 /3r(4/3) 
bj 


|36 J |V3r(4/3) 


V3i36 J | 1 /3r(4/3) 


^3j36 J | 1 /3r(4/3) 



5 





Pem 


1 







1 


{ 






5 


300 



FIG. 1. The induced emission probability P em as a function 
of the grid point number J for k„L = 10 (k/n n = 0.01, sech 2 
cavity mode profile). 




FIG. 4. The induced emission probability P em as a func- 
tion of the interaction length k„I for kj 'k„ = 0.01 (a) and 
k/K n = 0.1 (b) (J = 200, first excited sinusoidal cavity mode 
profile) . 





FIG. 2. The induced emission probability P em as a func- 
tion of the interaction length n n L for k/n n — 0.01 (a) and 
k/K n = 0.1 (b). The cavity mode is the sech 2 profile. Solid 
lines were obtained by analytical calculations, dot ones by our 
method (J = 200). 



FIG. 5. The induced emission probability P em as a function 
of the interaction length k„L for the sech 2 mode profile (a) 
and the Gaussian profile (b) (k/n n — 0.1, J = 300). 
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FIG. 3. The induced emission probability P em as a function 
of the interaction length n n L for k/K n — 0.01 (J = 100, 
fundamental sinusoidal cavity mode profile) . 
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